Dynamic facilitation picture of a higher-order glass singularity 
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We show that facilitated spin mixtures with a tunable facilitation reproduce, on a Bethe lattice, the 
simplest higher-order singularity scenario predicted by the mode-coupling theory (MCT) of liquid- 
glass transition. Depending on the facilitation strength, they yield either a discontinuous glass 
transition or a continuous one, with no underlying thermodynamic singularity. Similar results are 
obtained for facilitated spin models on a diluted Bethe lattice. The mechanism of dynamical arrest 
in these systems can be interpreted in terms of bootstrap and standard percolation and corresponds 
to a crossover from a compact to a fractal structure of the incipient spanning cluster of frozen spins. 
Theoretical and numerical simulation results are fully consistent with MCT predictions. 



Although the glassy state of matter has been a long 
time fascinating topic for physicists, it still presents sev- 
eral mysterious aspects Perhaps, the most controver- 
sial one is the very nature of the glass transition: that 
is the question of whether vitrification is a purely dy- 
namical process or rather the (dynamical) manifestation 
of a genuine thermodynamic amorphous phase. In spite 
of many progresses (for reviews, see [2461) the problem 
remains widely open. On one hand, there are notorious 
experimental and computational difficulties related to the 
exceedingly long equilibration times of macroscopic sam- 
ples, and to the possibility of detecting unambiguously 
the elusive 'amorphous order'. On the other hand, the- 
oretical modelling of glassy systems has made clear that 
slow relaxation processes are ubiquitous and may result 
from very distinct mechanisms. 

In a situation in which it is unknown how much inter- 
twined thermodynamics and dynamics are, it would be 
particularly advantageous from a methodological point 
of view to identify those peculiar glassy features that can 
be reproduced with no reference to specific energetic in- 
teractions. In this respect, facilitated spin models first 
introduced by Fredrickson and Andersen are particularly 
useful as they are constructed in such a way to have 
a manifestly uninteresting thermodynamics. Hence they 
allow to clearly disentangle, albeit in an arguably artifi- 
cial way, dynamic aspects from static ones. 

The first assumption of the dynamic facilitation ap- 
proach is that, on a suitable coarse-grained lengthscale, 
one can model the structure of a liquid by an assembly of 
high/low density mesoscopic cells which have no static in- 
teraction. Binary spin variables, taking on value ±1, can 
be simply assigned to these cells. The next crucial step is 
to postulate that there exists a timescale over which the 
effective microscopic dynamics takes a deceptively simple 
form: local changes in cells structure occur if and only 
if there is a sufficiently large number, say /, of nearby 
low-density cells (/ is called facilitation parameter). The 
latter assumption is actually difficult to derive by analyt- 
ical means. Nevertheless, it can be justified on physical 
grounds: it mimics the cage effect and gives arise to a 



variety of remarkable, and sometimes unexpected, glassy 
features, even if the thermodynamics is completely triv- 
ial [5(. Although this line of research has been pursued 
very actively in recent years, little attention has been 
paid to the possibility of reproducing, within this frame- 
work, more complex types of glassy behaviour (for some 
exceptions, see @, E|). In fact, even simple schematic 
MCT models predict the occurrence of topologically sta- 
ble singularities of higher complexity 2]. Liquids con- 
fined in a disordered porous matrix [10( and attractive 
colloids 111, |l2j are some examples in which these sce- 



narios have been recently observed. Solvable microscopic 
realizations of such systems would be highly valuable 
both for a deeper understanding of complex glassy fea- 
tures and for clarifying the limits of MCT. This is quite a 
delicate issue, however, as mean-field disordered systems 
which are supposed to be exactly described by MCT are 
spoiled by finite-size effects and pre-asymptotic correc- 
tions so strong to prevent a direct observation of MCT 



predictions [13|, [14| . For this reason, we follow here an 



alternative route. 

Motivated by recent findings on the annealed-quenched 
mixtures (Io| . we have generalized the dynamic facilita- 
tion approach by allowing for an inhomogeneous distri- 
bution of facilitation or, equivalently, as we shall see, of 
lattice connectivity. While the latter may originate from 
a geometrically disordered environment, e.g. a porous 
matrix, the former can be thought of as resulting from 
the coexistence of different lengthscales in the problem, 
e.g. mixtures of more or less mobile molecules/polymers 
with small and large size. By doing so, we show that 
facilitated spin models on Bethe lattice provide a close 
microscopic realization of the simplest higher-order bifur- 
cation singularity scenario predicted by MCT, in which 
two liquid-glass transition lines (of types A and B) join 
smoothly at a common endpoint. 

The model. Facilitated spin models consist of N non- 
interacting spins Oi = ±1, i = l,...,N with Hamilto- 
nian H = — fX^i *; evolving with a Metropolis-like 
dynamics: at each time step a randomly chosen spin 
is flipped with transition probability: w(o~i — > —o~i) — 
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min { 1, e _,l<Ti /' CBT } , if and only if at least / of its z neigh- 
boring spins are in the state —1 (hereafter h/k B = 1). 
On a Bethe lattice, the dynamics can be characterized 
by exploiting the relation with the bootstrap percola- 
tion 
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Il6j |. Disregarding the less interesting nonco- 
operative case, / = 1, in which there is no transition, 
one can distinguish two cases, i) For z — 1 > / > 1, 
the system undergoes a dynamical arrest: below a cer- 
tain temperature T c , the fraction of frozen spins $, which 
plays the role of the non-ergodicity parameter in MCT, 
jumps from zero to a finite value. This corresponds to 
the sudden emergence of a giant cluster of frozen spins 
(i.e. spins that are surrounded by more than / neigh- 
bouring 1 spins) with compact structure. This dynamic 
transition has a very peculiar hybrid nature: it is discon- 
tinuous and, at the same time, has diverging fluctuations 
as in continuous phase transitions. The geometric ori- 
gin of this behaviour has been understood quite in detail 
as being related to the divergence of the size of corona 



clusters near the transition 17, 181. Several results [1 



including those related to large scale cooperative rear- 
rangements responsible for slow dynamics [19( , have sug- 
gested a strong analogy with MCT. ii) For / = z — 1, z 
the transition is continuous as bootstrap percolation is 
equivalent to conventional percolation [151 ] ■ and for this 
reason it has attracted less interest. The facilitated spin 
dynamics of these systems has never been explored to our 
knowledge. To study the crossover between integer values 
of / we introduce a new facilitated spin model in which 
the facilitation strength can be continuously tuned. This 
is obtained by making the facilitation parameter a lattice 
site dependent quenched random variable, in close anal- 
ogy with the bootstrap percolation problem studied by 
Branco [2(|. For concreteness, we shall consider a Bethe 
lattice with coordination number z = 4, in which the 
facilitation /j is chosen from the probability distribution: 



P(fi) = (l-q) S fi<2 + (q-r) 5 fi<3 + r S ft 



(1) 



with 1 > q > r > 0. By tuning q and r we can 
thus explore the discontinuous/continuous glass transi- 
tion crossover. 

Exact results. The tree-like structure of the Bethe lat- 
tice allows for an exact calculation of the phase diagram. 
Following Refs. 15, 16, 2o| we get: 
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log (3q-l) 



, if < q < 1/2 
if 1/2 < q < 2/3. 



(2) 



The phase diagram is depicted in fig. [T] and comprises 
two lines that smoothly join at q = 1/2. As expected, the 
introduction of less facilitated spins increases the glass 
transition temperature, as compared to the pure case, 
q = 4j| ■ The precise nature of the two glass transitions 
depends on the behaviour of $ near T c . Denoting with 
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FIG. 1: Phase diagram for z — 4 and facilitation as in {T}. 
The dark region is the glassy phase. The dashed/solid line 
is the discontinuous/continuous transition. Inset: Fraction 
of frozen spins vs temperature for several values of q and 
r = 10" 3 . Below q = 1/2, $ jumps to a finite value at the 
transition which is represented by the dotted line. 



p = 1/(1 + c 1 / T ) the probability that a spin is 1 in 
thermal equilibrium, we find 
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$ = p[x 3 (4-3x)+6qx 2 (l-x) 2 + Arx(l - xf] (3) 
+ (l-p)[y 3 (4-3y) + 6qy 2 (l 

where y = p [x 3 + 3(1 — q)x 2 (l — a 
x = hp(l - 2q) + 
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y) 3 ] 

, and 



VV - 8p+ 12qp{l -p)j /2p(2-3g). 
The behaviour of $ as a function of the temperature is 
shown in the inset of Fig. [TJ for several values of q and 
with r = 10 -3 . Notice that $ generally depends on r 
while the phase diagram does not. For 0<<7<l/2,$ 
jumps to a finite value $ c = $(T C ) on the transition line, 
meaning that the infinite cluster of frozen spins has a 
compact structure. This structure turns out to be quite 
resilient against random inhomogeneities in the facilita- 
tion strength, in a rather large range of q. The critical ex- 
ponent /3 associated to the order parameter $ is obtained 
by expanding the above equations in the small parame- 
ter e = T — T c . As expected, we find $ — <i> c ~ e' 9 , with 
(3 = 1/2, that is the typical square- root dependence well 
known in MCT and in other systems with hybrid tran- 
sition. For 1/2 < q < 2/3 the glass transition changes 
nature as $ departs smoothly from zero with a power-law 
behaviour $ ~ e' 3 , and the correspondence with standard 
percolation suggests that the giant cluster of frozen spins 
has a fractal structure. Interestingly, in this case j3 de- 
pends on r: for r = one has (3 = 2, while as soon as a 
negligibly small amount of spins with /j = 4 is introduced 
into the system, r > 0, one has (3 = 1. The robustness 
of the latter behaviour reflects the fact that the mass 
of the fractal cluster of frozen spins is essentially dom- 
inated by the dangling ends, that is those parts of the 
cluster which are connected to the backbone by a single 
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FIG. 2: Persistence vs time for several temperatures, indi- 
cated in the key, and values of q. The transition is discon- 
tinuous/continuous for q smaller/larger than 1/2. We use 
r — 10 -3 , except for q — 0.52 where r = 0. The horizontal 
lines show the theoretical prediction for the plateau heights 
(the dashed one is the critical plateau at T = T c ). The sim- 
ulations are performed on lattices with N = 10 5 to 5 x 10 5 
sites and averaged over 2-20 samples. 



site. Only when r = 0, the latters are completely re- 
moved from the infinite cluster and the exponent changes 
to j3 = 2 [20(. Hence, the general scenario emerging for 
an arbitrary ternary mixture is that of two distinct glass 
transitions with (3 = 1/2 (for the discontinuous case), 
and j3 = 1 (for the continuous one). These critical ex- 
ponents reproduce exactly the MCT results for the F\2 
schematic model [2j, in which the memory kernel takes 
the form m\2 = Vi<p(t) + V2(f> 2 (t). 

Numerical simulations. We now turn to numerical sim- 
ulation to explore those features of equilibrium relaxation 
which are relevant for a comparison with MCT. Testing 
MCT in systems which are described by a somewhat ad 
hoc kinetic rule is particularly interesting because it al- 
lows us to probe the degree of universality of MCT results 
beyond the context (the actual Newtonian or Brownian 
liquid dynamics) and the approximations in which they 
were originally derived. The dynamics of facilitated spin 
systems is conveniently characterized by the persistence 
4>{t), i.e. the probability that a spin has never flipped 
between times and t. The long-time limit of <fi(t), 
which plays the role of the Edwards- Anderson param- 
eter in spin-glasses, is directly related to the fraction of 
frozen spins 4> = lim.t_j.oo 4>(t), which is known to describe 
broken ergodicity. Simulation results for the persistence 
are shown in fig. [2] for various q, above and below the 
glass transition temperature T c (q). To avoid dynamic re- 
ducibility problems the ternary mixture has a negligibly 
small fraction of spins with / = 4, typically r = 10~ 3 . We 




-2 2 4 6 



logio t 

FIG. 3: Equilibrium relaxation at criticality for discontinuous 
(q = 0.4) and continuous (q = 0.52) glass transition. Solid 
lines are power-law fits with exponent a = 0.224 and a = 
0.315. Inset: /3-relaxation time r vs temperature T. Solid 
lines are power-laws with exponent l/2a and 1/a, respectively. 



see that inside the glassy phase, the persistence attains 
a finite plateau which is in excellent agreement with the 
value of $ analytically computed in the previous section. 
On approaching T c (q) from above and for q < 1/2 we 
find typical signatures of MCT: two-step decay of 4>(i), 
late stage stretched exponential relaxation, and power- 
law form of equilibration time. When q — > 1/2 (and so 
— > 0), the range over which the time-temperature su- 
perposition principle holds shrinks, and consequently no 
simple scaling form of relaxation data at various temper- 
ature can be found in this limit. 

To make a more stringent test of MCT we now exam- 
ine in detail the critical dynamics on the glass transition 
lines, for which MCT predicts distinctive patterns of uni- 
versal critical behaviour. When the system is relaxing 
exactly at the critical temperature T c (q), MCT predicts 
the power law decay, 4>{t) — 4>(oo) ~ t~ a ( q \ irrespective 
of whether the long-time limit of persistence is finite or 
zero. Fig. [3] shows that indeed this power law behaviour 
is very well obeyed. A further remarkable prediction of 
MCT is that the exponent a(q) is intimately related, in a 
way that depends on the nature of the glass transition, to 
the power-law exponent describing the divergence of the 
characteristic time r on approaching the plateau at $ c 
(the so-called /3-relaxation regime): when temperature 
gets closer to T c (q) one has r ~ e -1 / 2 "^' for the discon- 
tinuous transition, and r ~ ^ 1 / a (i) for the continuous 
one. These two distinct behaviors are tested in the inset 
of Fig. [3] by assuming as a reasonable definition of the /3- 
relaxation time, r, the integral of <f>(t) from zero up to the 
time t* such that <fr(t*) = $ c . The estimated t behaves as 
a power-law of e near T c and, more importantly, the expo- 
nents perfectly agree with MCT, for both the discontinu- 
ous and continuous glass transition. The departure from 
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the plateau, the so-called a-relaxation regime, is also con- 
sistent with the exponent b obtained from the universal 
MCT relation T(l - af/T{\ - 2a) = T(l + b) 2 /T(l + 2b), 
though this regime is much more difficult to analyse due 
the exceedingly long time scales of numerical simulations. 

The excellent agreement we found with MCT could not 
be obviously anticipated by looking at the 'simple' facili- 
tation rule and is all the more remarkable by considering 
that i) no exact mode-coupling equation holds for facil- 
itated spin models 2l|] , and ii) the facilitation rule is 
admittedly quite remote from the actual microscopic liq- 
uid dynamics. Furthermore, our results imply that no 
thermodynamic transition is actually needed to observe 
a complex glassy scenario, and suggest that the difficulty 
to observe MCT predictions in mean-field spin-glasses, is 
generally related to the fact that the dynamic and static 
transition properties in these systems are tightly twisted. 

Conclusions. To summarize, we have extended the dy- 
namic facilitation approach to glassy systems by intro- 
ducing facilitated spin mixtures which exhibit a crossover 
from a discontinuous to a continuous glass transition, and 
shown that they provide a close microscopic realization 
of the simplest higher-order glass singularity predicted by 
MCT. Although we focused on Bethe lattice with fixed 
connectivity and random distribution of facilitation, it 
is worth to remark that one can consider an equivalent 
variant, in which the Bethe lattice is diluted and the fa- 
cilitation is uniform. For example, results qualitatively 
similar to those reported above can be obtained when 
the local lattice connectivity Zi is distributed according 
to P(zi) = (1 - q)5 ZiA + (q- r)5 Zi , 3 + r5 Zu2 , and fa = 2 
on every site. In this case, less connected sites turn out 
to have a smaller probability to flip, just as if they were 
less facilitated. We thus expect, and indeed find, that 
both variants give essentially the same results. Thus, the 
crossover between the two glass transitions is obtained 
by varying either the local connectivity or the facilitation 
strength, and corresponds to a passage from bootstrap to 
standard percolation transition. Consistently with MCT 
predictions, the critical dynamics on both glass transition 
lines is characterized by power-law decays with exponents 
closely related to those describing the divergence of the 
/3-relaxation time. We also mention that when a fraction 
of spins is overfacilitated as compared to the case /, = 2, 
i.e. when P(fi) = (1 - q)S fi:2 + (q-r)5f z: i+r6f t , , there 
is only a discontinuous glass transition line. 

Our work can be naturally developed in several di- 
rections. The most obvious one consists in considering 
kinetically constrained particle mixtures, what would be 
a much closer microscopic realization of a fluid confined 
in a disordered porous matrix. The comparison with ge- 
ometrically constrained lattice glass mixtures [22j , would 
then allow for a better understanding of the distinguish- 
ing features of glassy systems with and without an ideal 
Kauzmann transition. We then expect that introduc- 



ing a suitable static attraction [8| can account for non- 
monotonic dependence of relaxation time on attraction 
strength and, possibly, for phase reentrant behaviour. 
Finally, the existence of kinetic models with a glass tran- 
sition in d = 2 23], suggests that there is no apparent 
limitation to extending our approach to finite dimensions. 
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